Decoherence Free Subspaces for Quantum Computation 



D. A. Lidar/ a > I. L. ChuangW and K. B. Whaley( a ) 
(a) Department of Chemistry, The University of California, Berkeley, CA 94720 
(b) IBM Almaden Research Center, San Jose, CA 94120 

Decoherence in quantum computers is formulated within the Semigroup approach. The error 
generators are identified with the generators of a Lie algebra. This allows for a comprehensive 
description which includes as a special case the frequently assumed spin-boson model. A generic 
condition is presented for error-less quantum computation: decoherence-free subspaces are spanned 
by those states which are annihilated by all the generators. It is shown that these subspaces are 
stable to perturbations and moreover, that universal quantum computation is possible within them. 
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Decoherence remains the most important obstacle to 
the exploitation of the speedup 0] promised by quantum 
computers. To this end a remarkable theory of quan- 
tum error correction codes (QECC) has recently been 
constructed [0, in which a logical quantum bit (qubit) 
is encoded in the larger Hilbert space of several physical 
qubits. This "active" error- correction approach builds on 
the assumption that the most probable errors are those 
which occur independently to a few qubits during a rea- 
sonable time interval. However, correlated errors, which 
affect many or all qubits, may also be likely in some 
experimental realizations, particularly when qubits are 
physically close (for example, nuclear spins in a molecule) 
||. Such situations motivate the present study of an 
alternative "passive" error-prevention scheme, in which 
logical qubits are encoded within subspaces which do not 
decohere because of reasons of symmetry. The existence 
of such decoherence-free (DF) subspaces has been shown 
by projection onto the symmetric subspace of multiple 
copies of a quantum computer H , and by use of a group- 
theoretic argument || . Construction of these subspaces 
has been performed explicitly for certain collective error 
processes in the spin-boson model ™7|. In this Letter 
we formulate a general theory for decoherence in quan- 
tum computation (QC) within the powerful semigroup 
approach |^,^], and show that this provides a rigorous 
and comprehensive criterion for construction of DF sub- 
spaces for an arbitrary Hamiltonian. 
The Semigroup Approach — . The dynamics of a quan- 
tum system S coupled to a bath B (which together form 
a closed system) evolves unitarily under the Hamiltonian: 
VLsb = H ® I B + ts ® H B + H/, where H, H B and H/ 
are the system, bath and interaction Hamiltonians, re- 
spectively. I is the identity operator. In the Semigroup 
approach one shows that under the assumptions of (i) 
Markovian dynamics, (ii) "complete positivity" M, and 
(iii) initial decoupling between the system and bath JTo| , 
the following master equation provides the most general 
form for the evolution of the system density matrix p: 



a, /3=1 

L Fq ,f>] = [F a ,pFt] + [F a p,Ft] 



(2) 
(3) 



The commutator involving H is the ordinary, unitary, 
Heisenberg term. All the non-unitary, decohering dy- 
namics is accounted for by Lp. The time-independent 
Hermitian coefficient matrix A = {a a fj} contains the 
information about the physical decoherence parameters 
(lifetimes, longitudinal or transverse relaxation times, 
and various equilibrium parameters such as stationary 
polarization or magnetization) 

The {FqJ^q (F = I) constitute a basis for the vec- 
tor space of bounded operators acting on H, the N- 
dimensional system Hilbert space. This operator-space 
may be restricted - see the classification below. As such, 
the set {F a }*' =1 forms an M-dimensional Lie algebra £, 
with an N x N (generally M < N 2 — I) matrix repre- 
sentation {F a }^ =1 appearing in Eq. (Q) (we omit the hat 
symbol for matrices). Physically, the {F a }a=i describe 
the various decoherence processes: in the QC context 
they are the error generators. They are often determined 
implicitly by the interaction Hamiltonian: 



H, 



(4) 



dp 
dt 



[H,p]+L D [p] 



(1) 



where {B a } are bath operators (see Ref. |Tl[ for exam- 
ples) . 

Decoherence of a Quantum Register — . Consider a quan- 
tum computer made of K qubits. States in the corre- 
sponding TV = 2 K -dimensional register Hilbert space H 
are tensor products of single qubit states \e K ), e K = 0, 1. 
It is convenient to adopt the following classification of 
decoherence models of interest, in terms of the above Lie- 
algebraic scheme: (i) " Total decoherence:'''' This provides 
the maximum possible complexity of error generation, in 
which combined errors from any number of qubits are 
generated. As is well known, single-qubit errors can be 
fully described by the three Pauli matrices [i.e., the defin- 
ing representation of the Lie algebra su(2)]. Thus when 
\e K ) are the eigenstates of the cr% Pauli matrix, a sin- 
gle qubit can either undergo a phase-flip {cr z K ), a bit-flip 
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(cr^), or both (<r^). Taking into account also the pos- 
sibility of no single-qubit error, there are 4 possibilities 
per qubit, so that the maximal total number of com- 
bined errors on K qubits is M = A K — 1, if we disregard 
the case of zero overall errors. The Lie algebra su(N) 
has N 2 — 1 generators, so the corresponding M tensor 
products of Pauli matrices {F Q } form the defining rep- 
resentation of C = su(2 K ). (ii) 11 Independent qubit de- 
coherence:" In this, the ideal starting point for QECC, 
we have the much simpler case of merely one indepen- 
dent error per qubit, with all other qubits unaffected. 
There clearly are 3K such errors, each formed by taking 
the tensor product of a single Pauli matrix on one qubit 
with the identity on all the rest. Since errors on differ- 
ent qubits commute, this leads to a representation of the 
Lie algebra C — (B^ =1 su K (2). (iii) "Collective decoher- 
ence:" One could also consider the extreme case of all 
qubits undergoing the same decoherence process simul- 
taneously 0, i.e., assuming full permutation invariance 
of the qubits. There are then just 3 possible errors and 
C = su(2). (iv) 11 Cluster decoherence:" Situations inter- 
mediate between the above 3 cases follow when the reg- 
ister can be partitioned into clusters k of K' qubits, with 
collective decoherence taking place within each cluster, 
but the clusters decohering independently. This leads to 

C = ®kLi su k{2)- Lastly, a very interesting case (dealt 
with in detail below) arises when a symmetry (e.g., per- 
mutation invariance) is broken perturbatively. 
Conditions for Decoherence- free dynamics — . Within 
the extremes delineated by the above categorization, a 
particularly interesting question is: what are necessary 
and sufficient conditions for the existence of a generic 
DF subspace? By generic (as opposed to general), we 
mean that one should (a) avoid fine-tuning of the noise 
parameters characterizing the decoherence processes, and 
(b) avoid a dependence on initial conditions. Suppose 
that is a basis for an TVo-dimensional invariant 

DF subspace Iri C Ti. In this basis, we may express states 
as the density matrix 



To satisfy condition (b) above, each of the terms in paren- 
theses must vanish separately. This can only be achieved 
if there is just one projection operator \n)(m\ in each 
term. The least restrictive choice leading to this is: 
cf n — efSm. Eq. (JsJ) then becomes: 



l f q ,f« [p] 



N 

J2pij\i)(j\ ( 

ij=l 



/3* o 



„/3* a 



(7) 



Assuming cf ^ then yields: — + ^p- = 2. This has 

art 

to hold in particular for a = (3. With z = c a j/c a i, we 
then obtain z + 1/z* — 2, which has the unique solution 
2 = 1. This implies that c a i must be independent of i and 
therefore that F a \i) — c a \i), Va. As a result we conclude 
that: [FcF^jji) =0. If C is semisimple (has no Abelian 
invariant subalgebra) fll^| then the commutator can be 
expressed in terms of non-vanishing structure constants 

We 



fl p of the Lie algebra: [F a ,Fg] = 
then arrive at the condition on the structure constants 
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Va, (3. 



(8) 



Now, it is known that the structure constants themselves 
define the M-dimensional "adjoint" matrix representa- 



tion of C 12 



ad(F Q ) 



7.0 



p 

Ja,f3- 



Since the generators 

of the Lie algebra are linearly independent, so must be 
the matrices of the adjoint representation. One can read- 
ily show that this is inconsistent with Eq. (@) unless all 
c 7 = 0. We have thus proved Jl3| : 

Theorem 1 A necessary and sufficient condition for 
generic decoherence- free dynamics (Ld[/3] = 0) in a sub- 
space Ti. = Span[{|i)}^ > 1 ] of the register Hilbert space, is 
that all basis states \i) are degenerate eigenstates of all 
the error generators {F Q }: F Q |i) = c a \i), Va; or, if C is 
semisimple, that all \i) are annihilated by all {F Q }: 



No 

p = E ^'Woi 



(5) 



Va, i. 



Consider the action of the error generators on the basis 
states: F a \i) = Ylj=i c ?j\j)- The DF dynamics condition 
is Ld[p] = 0, so that by Eq. (|l|) the dynamics is purely 
unitary in the subspace H. Consider then Eq. (^): con- 
dition (a) above implies that each of the terms Lf ctj f 3 [p] 
should vanish separately Va,f3. A straightforward calcu- 
lation yields: 



l f q ,fJ/5] 
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ij,mn=l 



Pij (2c£,cg,|n)<m| 

CHOI - c fm C «mN}H) 



(G) 



(9) 



Equivalently, the DF subspace is spanned by those states 
transforming according to the 1-dimensional irreducible 
representations (irreps) of the Lie group with algebra C. 
Those states are singlets. The size of the DF code pro- 
vided by this subspace is its dimension No, which can be 
used to further encode log 2 (A r o) logical qubits. 

Note also that by Eq. (g): Hj|i) <g> \b) = 0, where |6) 
is any bath state. Theorem 1 thus not only reduces the 
identification of DF subspaces to a standard problem in 
representation theory of Lie algebras, but also has the ex- 
pected physical interpretation, namely that the DF states 
are those that are annihilated by the interaction Hamil- 
tonian. (Note that this is only a necessary condition.) 
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Effect of the system Hamiltonian — . While p, by con- 
struction, is unaffected by the error generators, the ab- 
sence of decoherence may still be spoiled by the system 
Hamiltonian itself. To see this explicitly, consider the 
mixed-state fidelity: 

F(t) = Tr[p(0)p(t)] = Tr [p(0) exp(L t) [p(0)]} , (10) 

which is a natural measure of the decay of quantum co- 
herence due to coupling of the system with the environ- 
ment. In ideal quantum computation, one would like to 
have F(t) = 1, corresponding to perfect, noiseless mem- 
ory. In reality F(t) = 1 — e, e > 0. A formal power 
expansion yields: 

oo oo 1 / , x n 

F(t) =£-Tr[p(0)(Ln P (0)]]^- - , (11) 



n=0 



n=0 



where the "decoherence times" are: r n = 
{Tr[p(0) (L)" [p(0)]]}~ 1/n . In particular the first order 
decoherence rate is: 



l = Tr[p(0) Lp(0)] 



(12) 



Since Tr p[H, p] = (by cyclic permutation), it thus 

follows from Eq. (|l|) that 1/ri = for p. However, as is 
easily checked, generally 1/t2 ^ because H may cause 
transitions outside of TL. Therefore, the full dynamics in 
H, including the effect of the system Hamiltonian, is DF 
to first order. 

Effect of Symmetry Breaking Perturbations — . Suppose 
we have identified the DF subspace for the Lie algebra C 
underlying hry. Let us consider the effect of adding new 
error generators {G p }^ =1 which perturbatively break the 
symmetry, i.e., which do not belong to C. We assume 
that the {G p } are due to an additional interaction Hamil- 
tonian Hj which can be identified as appearing with a 
small parameter e in the full system-bath Hamiltonian: 
flsB = H + Hb + Hj + ell'j. Then the new terms added 
to Lu are: 

M P 

L 'dIp\ = ( a "P L F Q , e G p [p] + a* ap L tGp ,F a [p\) 

a—1 p—1 
P 

+ a P? L «G p ,eGjp]. (13) 

p,q=l 

Under the assumption e <C 1 we may neglect the last 
term since it is 0(e 2 ). As for the terms in the double 
sum, F Q p = pFjj = by Eqs. @ and (||). Expanding 
out the remaining terms leaves: 



M 



*>d\p\ ~ ^Yl ««,pP G p F « + H.c. 



(14) 



While this will generally take the singlet states outside 
of the DF subspace, this effect is also readily seen to be 
only of second order, because the first-order decoherence 
time [Eq. (p~2f)] is now given by: 



i p 

7 =^{« a , p Tr[^(0)p(0)GtF„] 
Tl P =i 

+ < p Tr[p(0)FtG p p(0)]}=0, 



(15) 



by cyclic permutation under the first trace. The higher 
order decoherence-times, r„, clearly involve e™ and can 
thus be made negligible. Therefore we have proved that 
the DF subspace is stable to first order under a 
symmetry breaking perturbation. 
This property is very promising from a quantum com- 
putational perspective, since one should be able to ap- 
ply standard QECC techniques to correct errors which 
then occur within the DF subspace. Of particular con- 
cern are errors which take states out of the DF subspace; 
these are analogs of amplitude damping errors, which ab- 
stractly model for example scattering and spontaneous 
emission processes. Such errors can be corrected by sim- 
ple codes Q, for example, by taking the DF singlet 
states as the computational basis states, and combin- 
ing them into QEC codewords. Provided that causes 
independent errors on different singlet states, we can con- 
clude from the threshold theorem (ljjl^l that as long as 
e is sufficiently small, the QECC encoding will render 
quantum computation within Tt robust against these er- 
rors. Typical estimates of the threshold error probability 
range from 10~ 6 to 1CT 3 Jl6| and are extremely difficult 
to achieve in practice. The error probability is usually 
proportional to e 2 . However, within Tt!, the error proba- 
bility is reduced to e 4 . Thus, QC within a DF subspace 
has potentially significant advantages. 
The dimension of DF subspaces: the size of codes — . As 
shown in Rcf. |tJ for the spin-boson model, in the limit 
of collective decoherence (i.e., when C = su(2)) the size 
of the DF subspace is: 



N K -^K-hog 2 K. 



(16) 



a—l p—1 



The encoding efficiency Nq/K, is thus asymptotically 
unity. However, in the opposite limit of independent 
qubit decoherence, C = ©^L 1 su K (2), which is addressed 
by QECC, there does not exist a DF subspace |oL The 
size of the code obtained in intermediate cases of cluster 
decoherence can be estimated from Eq. ( |l6| ) by replacing 
K by K' (the number of qubits per cluster), as long as 
K' < K . However, the most interesting situation arises 
in the perturbative scenario. Imagine a case of collective 
decoherence symmetry which is perturbatively broken by 
small independent couplings between individual qubits 
and the bath. As long as the symmetry-breaking inho- 
mogeneities are not too strong, we can conclude that, to 
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first order, the exponentially large DF subspace is still 
available. 

Universal Quantum Computation — . Our discussion so 
far has centered on the preservation of quantum mem- 
ory. To complete it we still need to show that uni- 
versal quantum computation can actually be performed 
in the DF subspace. As is well known, the controlled- 
not operation, together with arbitrary single qubit ro- 
tations, can generate any unitary operation [p|. The 
corresponding unitary operations are implemented by a 
driving Hamiltonian H^, which contains experimentally 
manipulable, time varying parameters, together with the 
system Hamiltonian H. 

We now give an example of universal 1- and 2-qubit op- 
erators acting on a 4-dimensional singlet subspace. Let 
< i < 3 be singlet states. These 4 states span 
2 encoded qubits | <?i ) , | <?2 ) where qiq2 (with qj = 0,1) 
is the binary representation of i. A controlled-NOT 
gate can be constructed from a Hamiltonian represented 
in the encoded basis by the following combination of 
projection operators: H^ ot = c{t) [|11)(10| + 1 10) <1 1 1] . 
Here c(t) is a time-dependent classical control pa- 
rameter. Upon exponentiation this yields the fa- 
miliar conditional unitary operator form. Single 
encoded-qubit rotations can be constructed from, e.g., 
H™* = no(t)[|Oi}(00| + |00)(01| + |ll)(10| + |lO)(ll|]. 
The generalization to larger singlet space systems is 
straightforward: one constructs the appropriate projec- 
tion operators on the singlet states. By construction 
the resulting gates will leave the dynamics DF. Thus, in 
principle, universal DF-QC is possible within the singlet 
subspace. The main experimental challenge will involve 
implementation of the corresponding operations on the 
physical qubits. In addition, one should expect the actual 
implementation to involve some of the amplitude damp- 
ing errors discussed above, i.e., some H<j operations will 
take the singlets out of the DF subspace. However, as 
long as QECC is invoked, our previous arguments show 
that DF-QC is still possible. 

Conclusions — . It was shown how decoherence in QC 
can be described very generally in terms of the Semigroup 
approach. The usual QC "error generators" were identi- 
fied with the generators of a Lie algebra, whose identity 
depends on the pertinent decoherence process. Without 
reference to a specific system-bath interaction model, we 
derived the generic condition for DF subspaces: these are 
spanned by those states which are annihilated by all the 
error generators. We showed further that the DF sub- 
spaces are stable to first order under symmetry breaking 
perturbations, which allowed us to extend their utility 
by application of QECC. Finally, we showed that the DF 
subspaces support universal quantum computation. 
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